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ADAM NOWAK AND KRZYSZTOF STEMPAK 

Abstract. We establish a symmetrization procedure in a context of general orthogonal 
expansions associated with a second order differential operator L, a 'Laplacian'. Com- 
bined with a unified conjugacy scheme furnished in our earlier article it allows, via a 
suitable embedding, to associate a differential-difference 'Laplacian' L with the initially 
given orthogonal system of eigenfunctions of L, so that the resulting extended conjugacy 
scheme has the natural classical shape. This means, in particular, that the related 'partial 
derivatives' decomposing L are skew-symmetric in an appropriate L 2 space and they com- 
mute with Riesz transforms and conjugate Poisson integrals. The results shed also some 
new light on the question of defining higher order Riesz transforms for general orthogonal 
\ expansions. 

-I— > 

1. Introduction 

The seminal article of Muckenhoupt and E. M. Stein [2] initiated the investigation 
of conjugacy for discrete and continuous nontrigonometric orthogonal expansions. In the 
recent years a considerable activity could be observed in studying conjugacy, or better Riesz 
transforms, for orthogonal expansions in one-dimensional and multi-dimensional settings 
related to general second order differential operators. 

A variety of papers has been devoted to the study of objects being ingredients of con- 
jugacy notions defined by different authors in many particular situations. In connection 
to a dynamic development of investigation of conjugacy problem in different settings, a 
natural demand appeared on a general and universal definition of Riesz transforms. The 
authors' paper [3] was an attempt to provide a reasonable answer to this important de- 
mand by offering a unified conjugacy scheme that includes definitions of Riesz transforms 
and conjugate Poisson integrals for a broad class of expansions. The postulated definitions 
were supported by a "good" L 2 -theory, existence of Cauchy-Riemann type equations, and 
numerous examples existing in the literature which are covered by the scheme. 

There is, however, a shortcoming of this unified conjugacy scheme manifested in a lack 
of symmetry in the decomposition 

d 

3=1 

of a given second order partial differential operator L, a 'Laplacian', acting on functions 
on a d- dimensional domain X = (b, c) d , oo < b < c < oo. Here A > is a constant, 
5j are 'partial derivatives' associated to L (first order partial differential operators, tr- 
acts on the jth coordinate), and 5* are their formal adjoints in an appropriate L 2 sense. 
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Riesz transforms of first order defined in [3] are formally given by Rj = SjL- 1 / 2 (or by 
Rj = 5jL -1 ' 2 IIq, see [3] for details), but a replacement of Sj by 5* in this definition is, in 
general, inappropriate since it may result in an operator taking L 2 functions out of L 2 . 

Asymmetry of the decomposition of L has, in fact, a deep impact onto the whole con- 
jugacy scheme postulated in [3]. To be precise, taking into account existing examples, it 
seems that the case of an operator L acting on the whole space M. d is not really affected 
by this asymmetry Therefore, in what follows we consider only the case (b, c) ^ M, and 
assume, without any loss of generality, that = b < c < oo. Then a possible way of 
overcoming the lack of symmetry is provided by a symmetrization procedure, which is the 
purpose and the main achievement of the paper. This procedure is to some extent inspired 
by a situation of certain orthogonal systems appearing in the theory of Dunkl operators, see 
Section [3] for more comments. 'Partial derivatives' emerging from the symmetrization pro- 
cedure, contrary to Sj's, are skew-symmetric as it happens in many classical cases including, 
in particular, the usual Euclidean partial derivatives, Dunkl operators, left-invariant vector 
fields on Lie groups, etc. 

Throughout the paper we use a fairly standard notation. The symbols A and dx will 
always refer to the Euclidean Laplacian, A = £) i=1 °\ , and Lebesgue measure acting, or 
considered, on an appropriate domain in M. d like, for instance, = (0, oo) d . The symbol 
N is used to denote the set of nonnegative integers, N = {0, 1, 2, . . .}. Finally, (•, denotes 
the canonical inner product in an appropriate L 2 space, where \x is a given measure. 



2. Initial situation 

Our starting point is the situation discussed in [31 Section 2], where the concept of 
studying conjugacy for orthogonal expansions is based on the existence of a second order 
differential operator playing a similar role to that of the standard Laplacian in the classical 
harmonic analysis. Below, d > 1 will always denote the dimension. 

We first consider the following one-dimensional objects, which in a while will serve as 
building blocks of <i-dimensional product structure: 

• an open (possibly unbounded) interval IcK; 

• a system {//j : i = 1, . . . , d} of absolutely continuous measures on X, ^{dxi) = 
Wi(xi)dxi with strictly positive densities Wi G C 2 (X); 

• a system {Li : % — 1, . . . , d} of second order differential operators defined on C 2 (X). 

Here, in this paper, we exclude the case X = M (see the comment in Section [TJ, so without 
any loss of generality we may assume that X = (0, c), < c < oo. For each of the operators 
Li, in order to ensure existence of the associated 'derivative', we assume the decomposition 

(2.1) L 4 = a, + ^, 

where is a nonnegative constant, and Si is a first order differential operator (a 'derivative') 
of the form 

d 

Si = Pi(xi)- h qi(xi) 

OXi 

with real coefficients Pi G C 2 (X), q { G C l (X), Pi(xj) ^ for G X; here S* represents 
the formal adjoint of Si in L 2 (X,Hi), 

w f -(xi*) 
$i = -Pi(xi)-z- + Qi( x i) ~ Pi( x i) h \ ~ Pi( x i) 

OXi Wi(Xi) 
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determined by the identity 

Notice that 5* 7^ —5. Thus, a posteriori, each Lj is a linear operator with continuous 
real- valued coefficients and negative leading term coefficient, 

d 



t 2/ ^ 9 



2 Pl (xM(^+ti(^ ir:iJ '' ] 



Wi(xi) \ dx 



+ Qi( x i) - (Pi{xi)qi{xi))' - Pi(xi)qi{xi) W% X \ + a { . 

Wi{Xi) 

Moreover, (12. ip implies that each Lj is symmetric and nonnegative on C 2 (X) C L 2 (X, fa). 

Now we are in a position to specify a <i-dimensional setting that is suitable for further 
development. We equip the space X = X x . . . x X (d times) with the product measure 

// = fa <g> . . . ® fa. 

We consider the <i-dimensional 'Laplacian' L (more precisely, L is a generalization of —A) 
defined initially on C 2 {X) by 

L = L l + . . . + L d , 

where each Li is understood as a one-dimensional operator acting on the ith axis. Note 
that in view of the previous assumptions, L admits the decomposition 

d d 

L = A + <^> where A = > 0; 

4=1 4=1 

here the indices of 5 and 5* indicate also on which axes actions of these operators take 
place. 

Next, we introduce an orthogonal system associated with L. With no loss of generality 
we may restrict to L 2 -normalized systems. Assume that for each i = l,...,d, there 
exists an orthonormal and complete in L 2 (X,fa) system {(p k ^ : kj G N} consisting of 
eigenfunctions of Li, with the corresponding eigenvalues {A^? : h G N}, i.e. Li^ k . = 
Vk- ■ Here and below we assume for simplicity that <p^, G C°°(X), but in fact much less 
regularity is needed (we omit a discussion in this direction since it could affect the main 
line of thought of the paper). For a multi-index k = (k\, . . . , kd) G N d we define 

(1) o o (d) 

Then {ip k : k G N d } is an orthonormal basis in L 2 (X,fa) consisting of eigenfunctions of L, 
Lip k = X k ip kl where X k = \£ + . . . + A^. 

In addition, ip k G C°°(A'). 

We impose the following technical assumptions on the systems of eigenvalues and eigen- 
functions, which seem to be unavoidable on the considered level of generality. For every 
i — 1, . . . , d, we assume that the one-dimensional eigenvalues are indexed in the (strictly) 
ascending order, Aq < A^ < \% < ■ ■ ■, and lim^ X k = 00. Consequently, the set 
{X k '■ k G N d } of mult i- dimensional eigenvalues may be arranged into an increasing and 
divergent sequence 

A < Ai < A 2 < . . . , A m — > 00 when m — > 00. 
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Moreover, we require the associated 'partial derivatives' Si to be Inconsistent with the 
orthogonal system, i.e. for each i 

and 

(^gWgU = mvll^)^ h,m g N. 

All these assumptions are not too restrictive, as may be seen by various examples given in 
[3j Section 7]. 

In the situation described above it is not hard to check that the 'Laplacian' L is sym- 
metric and nonnegative on C 2 (X) C L 2 (X,fi), and the constant A in the decomposition 
of L does not exceed the smallest eigenvalue, A < A , cf. [3J Lemma 1]. However, from 
the conjugacy point of view, the following fact is essential (see [3J Lemma 2]): given 
i = 1, . . . ,d, the 'differentiated' system {Siipk '■ k G N d } is orthogonal in L 2 (X, /z); further- 
more, ll^i^jkH^a/^ ^ = — dj. This fact leads to investigating also 'Laplacians' standing 
behind the systems {5i(pk}, i — 1, ■ ■ ■ ,d, and this turns out to be a crucial point in con- 
structing proper conjugacy scheme for general orthogonal expansions. Define 

Mj = A + SjS* +J2 S *i S i = L + MIL j = l,...,d, 
where [5^, <5*] is the commutator 

[Sj, 6*} = 5j5* - 5*5j = 2 Pj {x j )q' j {x j ) - Pj {xj) 

By the very definition it follows that each Mj is symmetric and nonnegative on C 2 (X) C 
L 2 (X,fi). Moreover, for each j = l,...,d, the system {8jipk : k G N d } is an orthogonal 
system of eigenfunctions of Mj, with the corresponding eigenvalues {A^ : k G N d }, i.e. 

Mj(8jipk) = Xk(Sj^k); 

see [21 Lemma 5]. The operators Mj (or rather their suitable self-adjoint extensions) are 
used to generate so-called modified Poisson semigroups that play an important role in the 
conjugacy scheme for orthogonal expansions proposed in [3], see [31 Section 5] for details. 

The main inconvenience of the theory postulated in [3J is a lack of symmetry in principal 
objects and relations, and this phenomenon has roots in the asymmetry between the 
'derivatives' 5j and their adjoints 5*. In consequence, definitions and the conjugacy scheme 
established in [3J admit essential deviations from the classical shape, see [3j Sections 5,6]. 

The main idea of this paper is to overcome the problem by embedding the situation 
considered in [3] into a more general setting, where the associated derivatives are skew- 
symmetric and the related conjugacy scheme has precisely the classical shape. The price is, 
however, that the related extended 'Laplacian' and 'derivatives' are differential-difference 
operators rather than purely differential ones. It is remarkable that most definitions and 
relations in the setting of [3j may be then recovered by suitable 'projecting' from the 
extended situation. However, in some cases the projection procedure leads to different 
and seemingly even more natural definitions. This remark concerns especially higher order 
Riesz transforms. 



Wj[Xj) v 



Wj(Xj) 
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3. Symmetrization 

We now describe the symmetrization procedure and the resulting symmetrized situ- 
ation. The construction is motivated to some extent by the setting of the Dunkl har- 
monic oscillator with the underlying reflection group isomorphic to 7L\ = {0, l} d ; we refer 
to [1] for more details concerning the Dunkl setting. Recall that X = (0, c) for some 
< c < oo and X = X d , and consider the space X = -Xsym x . . . x Xsym (d-times), 
where Xsym = (— c, 0) U (0, c). Notice that X is isomorphic to each of the 'Weil chambers' 
generated in X by reflections perpendicular to coordinate axes. We extend the measure \i 
to X by even extension of the one-dimensional densities W{, Wi(-Xi) = Wi(xi), Xi > 0; we 
keep using the same symbols for the extended objects. Further, we extend the coefficients 
of L by letting 

Pi(-Zi) = Pi(xi), Qi(-Xi) = -qi(xi), Xi > 0; 

again the emerging extended objects, including L, Mj and Sj defined by means of the 
extended coefficients, are denoted by still the same symbols. 

Definition 3.1. For a suitable function f on X define its 'partial derivatives' 
Djf{x) = pjixjj — ix) + qjixj) J — 



+ 



Pj i X i ) -77-T + Pj ( X J)- % \ X 3 , 

w j\ x j) 



f(x) - f((TjX) 



where <jj denotes the reflection in X in the hyperplane orthogonal to the jth coordinate 
axis, o~j{x\^ . . . , Xj, . . . , X(i) ("^l; • • • j •Ejj • • • j 

The result below follows by integration by parts and some elementary manipulations. 

Proposition 3.1. The operators Dj, j = 1, . . . , d, are skew-symmetric in L 2 (X, /i), D* = 
—Dj, in the sense that 

(D j f,g), = -(f,D J g)„ f,geCl(X). 

This motivates the definition of the extended 'Laplacian' L as 

d 

(3.1) h = A-J2Dl 

i=i 

Then each Di commutes with L, which is an important feature at this point. 

To state the next result it is convenient to introduce the following terminology. Given 
e e Zj, we say that a function / on X is e-symmetric if / o a,j = (—l) Sj f, j — 1, . . . , d. If 
/ is e-symmetric and Ej = (e 3 - = 1) then / is said to be even (odd) with respect to the 
joth coordinate. 

Proposition 3.2. The operator L is symmetric and nonnegative on C^(X) C L 2 (X, /i). 
Moreover, for any e-symmetric function f G C 2 (X) ; e G Z^, we have 

::,/' Vf ■ E W+ E W 

U*j=0} {j-ej=l} 

In particular, L/ = Lf when f is even with respect to all coordinates, and L/ = Mjf if f 
is odd with respect to the jth coordinate and even with respect to the remaining coordinates. 
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Proof. The first part follows from the decomposition of L in terms of the Dj. Justifying 
the remaining part may be done by computing the explicit form of L, which is 



hf(x) = Af(x)-J2\pKxi) 



+ ( qi(xi) 



Pi(X 



,97 
dx\ 

w[{xi) 



x) + 



2p i (x i )p' i (x i ) +p 2 i(xi 



Wi(Xi) 



x 



+ Pi{Xi 



Pi[Xi 



Wi(Xi) 
w'^Xi) 



+ Pi( x i) 



Qi( x i) )f(x) +p i (x i )q' i (x, 



df_ 

dxi 

f(x) + f((TiX) 



wAXi 



+ Pi( x i) ~ Qi( x i 



' f(x) - f(a l x) 



The proof is finished by comparing the above expression with the explicit forms of 5*5i and 
5i5*, which may be read off from the explicit expressions for L and Mj, see Section □ 



Next we extend the eigenfunctions (p k to X by letting ip^(-Xi) = ip k l '(xi), Xi > 0, 
i = 1, . . . , d. Then, automatically, given e G Zg, the function 8?... 5 £ d d (fik is e-symmetric 
by the way of extending the coefficients of 8j. It turns out that these are eigenfunctions 
of L. 

Lemma 3.3. Let e G 7L\ be fixed. Then 

L(5f . . . 5 £ / Vk ) = \ k (5? . . . 5 e /<f k ), k g N d . 

Proof. Combine Proposition 13.21 with the product structure of Sf 1 . . . 5 £ d d ip k and the fact 
that in the one-dimensional setting ip k is an eigenfunction of L = a + 5*5 and 5tpk is an 
eigenfunction of M = a + 55* (to be precise, the last fact was already invoked from [3] in 
the non-extended setting, but it easily carries over to the extended situation by the way 
of extending the coefficients of 5 and 5*). □ 

Note that for a given s G it may happen that for some k G M d the function 5 ei . . . 5 6 d d ipk 
vanishes identically. This occurs precisely when there is an % G {1, . . . , d} such that £j = 1 
and di = Xq , and k G is such that hi = 0. 

To construct an orthonormal system {$„} associated with L and related to the original 
system {(p k }, it is natural to consider first the one-dimensional case. Then the relevant 
multi-dimensional system will be obtained simply by taking tensor products. The con- 
struction below is partially motivated by the case of the classical trigonometric system. 
Let 



,(0, 



,(0 



rii even, 



1/2 5j i] (x-) 



v/2V>i+l)/2 "0 

In this place it seems to be natural to require that for each i 

,0 



n, 



odd. 



1, . . . , d, the derivative 

5i annihilates the first eigenfunction, S^q' = 0. This is equivalent to assuming that the 
constant a, from the decomposition of Lj is equal to the first eigenvalue, a, = Aq . In 
the multi-dimensional setting the requirement means the equality A = A . We emphasize 
that this is indeed the case, up to a convention explained in a moment, of all the classical 
examples given in [21 Section 7]. Consequently, does not enter the definition of 

above. On the other hand, notice that all the $n,- are well-defined and non-vanishing. 
By the facts mentioned earlier (cf. [31 Lemma 2]) each of the systems {$>n} : rij G N}, 
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i — 1, . . . , d, is orthonormal in L 2 (X SY m, For a multi- index n = (n 1; . . . , n d ) G N d we 
define 

$„ = $«®...®$g. 

The multi- dimensional system {$„ : n G N d } is orthonormal in L 2 (X, /i). Moreover, the 
$ n are eigenfunctions of L, as stated below. 

Lemma 3.4. We have 

L$„, = + . . . + Aj^J^n, n G N d , 

where |_-J denotes the integer part function (the floor function). 

Proof. Given e G Zg, notice that, up to a constant factor, $2fc- £ coincides with Sf 1 . . . 6^ d (fk 
whenever 2k — e G N d . Then Lemma [3.31 gives the desired conclusion. □ 

In what follows, for multi-indices n G N d we will use the notation 



(n) 



ni + 1 



^ + 1 



and (in particular) (n) = L^^J when n is a number. Then we may write shortly 

L$ n = A (n> $ n , n G N d . 

The 'real' picture that emerges from the above procedure may be then turned into a 
'complex' one. Indeed, define first in dimension one 

*g = 7f + is S nn * ^tl-i)> n* e Z, 

and then for a multi- index n = (ni, . . . , n^) G Z d , 

An easy argument shows that the system {\l/ n : n G Z d } is orthonormal in L 2 (X, //) and 
consists of eigenfunctions of L, 

L^ n = A H ^ n , n G Z d , 

where |n| = (|ni|, . . . , [ri^ j ) . 

We remark that the choice of signs in the construction of {$ ra } is in principle arbitrary. 
Our particular choice is motivated by the fundamental example below. 

Example 1. The basic example here (and in some sense a prototype) is the case of 
classical trigonometric expansions. Let d — 1 and consider the interval X = X = (0, it) 
equipped with the measure fi(dx) = -dx. Further, consider the one-dimensional standard 

Laplacian L = — -j^ on (0, it) and the related orthonormal basis in L 2 (X,n) of cosines, 



k = 0, 
cos/cx, k > 0. 



Clearly, Ltp k = Xkfk, where A^ = k 2 . In addition M = L = — 4^. Applying the sym- 
metrization procedure we arrive at the trigonometric systems 

{$ n : n G N} = | -j=, sin x, cos x, sin 2x, cos 2x, . . . j 
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and 

{ty n : n E Z} = j— ^=exp(inx) : n E zj 

on the interval X = (— 7T,7r). These are orthonormal bases in L 2 ((— 7r,7r), ^dx) of eigen- 
functions of the Laplacian L = — considered on (— 7T, 7r), L$ n = (n) 2 $ n , L\l/ n = n 2 \l/ n . 
This example may be easily generalized to arbitrary dimension d > 1. 

In the situation of Example 1 we could as well choose as the initial system {(pk} the 
system of sines. This, however, leads to a small obstacle since then the constant in the 
decomposition of L does not coincide with the first eigenvalue, as required above. On the 
other hand, the system of sines is commonly enumerated by k — 1, 2, . . ., excluding k = 0. 
To overcome these problems, we introduce the following technical convention: in the case 
just described, and also in similar cases as those of Fourier-Bessel systems (see [H Section 
7.8]), we formally treat A = and ip = as the first eigenvalue and the corresponding 
eigenfunction, respectively. Then we are in a position to apply the symmetrization, which 
leads to an extended system {<3? n } with $ = to be neglected (thus, in fact, is 
enumerated by n = 1, 2, . . ., as is the initial system). Clearly, the convention just described 
in dimension one induces an analogous convention in the multi-dimensional situation. 

Applying this convention to (pk(x) = \/2 sin kx, k — 1, 2, . . ., and passing to symmetriza- 
tion we arrive at the trigonometric system {& n } as in Example 1, but with $ — Vv2 
excluded. Notice that this time the extended system is not complete. This indicates that 
the symmetrization applied to the system of cosines provides a more natural way of em- 
bedding the system of sines into the extended symmetric situation. In other words, it is 
more natural to view the sines as the 'differentiated' system rather than the initial one. 



4. RlESZ TRANSFORMS AND CONJUGACY 

In this section we investigate the symmetrized setting from the conjugacy point of view. 
We define Riesz transforms and conjugate Poisson integrals associated to the extended 
'Laplacian' L, and then show that these definitions fit into a consistent conjugacy scheme. 
This scheme, including Cauchy-Riemann type equations, has precisely the classical shape. 
When the convention described at the end of Section [3] is in force, then the results below 
should be understood accordingly. 

Following a general concept, we define formally the Riesz transforms of order N > 1 by 
R l = D l £~\ 1 " 2 , \l\ = N; all necessary notions will be explained momentarily. To make this 
definition strict, we need to specify a suitable self-adjoint extension of L. Consider the 
operator 

(4-1) £/= £A (n> 

neN d 

defined on the domain 

(4.2) Dom£={/eL 2 (X,/i): £ \\ {n) (f, $ n ) M | 2 < 00}. 

We denote by M = (span{$ n : n E ~N d , A( n > 7^ 0}) 1 - the null subspace of £. Note that 
M is not necessarily trivial. Independently of the case, we shall always write IIo for the 
orthogonal projection of L 2 (X, /i) onto Af x . 
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Lemma 4.1. The inclusion C 2 (X) C Dom£ holds, so that £ is a nonnegative self-adjoint 
extension of the operator n L defined initially on C 2 (X). Moreover, the spectrum of £ 
satisfies 

{A , Ai, . . .} C a{&) C {0} U {A , A 1; . . .}. 

Proof. Here arguments are similar to those from the proofs of [3j Lemma 3] and [3l Lemma 
6]. We omit the details. □ 

The spectral decomposition of £ may be written as 

oo 

£/ = ^A m P m /, / G Dom£, 

m=0 

where the spectral projections are 

V m f= (f,®n)^n, 171 EN. 

{neN d :A (n) =A m } 

Next we define more strictly the Riesz transforms of order N > 1 by 
R i = D i2-m u ^ l = ; k) e N ^ {(0) } 0)}> 

where \l\ = l\ + . . . + Id = N is the order of the transform, and D l = D 1 ^ . . . D 1 ^ (since 
the Dj commute, any composition of them may be written in such a form). Notice that 
for the order one, if I = Cj (the jth coordinate vector), then D l = Dj and consequently, 
R l coincides with ZJj£~ 1//2 Ho; in what follows we will denote these operators by Rj, j = 
1, . . . , d. If N > 2 and |/| = N, it is customary to call the operators R l the Riesz transforms 
of higher order. To provide a fully rigorous definition of R l we use the spectral series of £ 
and set 

(4.3) R l f= {\n)Y m {f^n),D^ n , /GL 2 (X,/i); 

(notice that A/ n \ = may happen only when n = (0, . . . , 0)). To show that this formula 
indeed gives rise to L 2 -bounded operators we first need to have a closer look at the action of 
the 'derivatives' on the eigenf unctions. Recall that = dj, j — 1, . . . , d, and so Ao = A. 

Lemma 4.2. Given j — 1, . . . , d and N > 1, we have 

r)Nfo = H- 1 ) N/2 (rfh- a i) N/2 *n> Neven, 

with the convention that $ n = ifn(£N d . 

Proof. Because of the product structure we may and do assume that d = 1 (thus k and n 
are nonnegative integers). Recall that 

®2k = -TEVk, $2fc-l = 7= 7= Otp k . 

Since tfk = a/2$2A: is an even function we have 
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and since — D 2 = L — a and (L — ajipk = (A^ — a)(pi» (see Lemma f3 . 3 j) . 

D$ 2Jfc _i = —D 2 ( j=L $ 2fc ) = (L - a) ( j=t $ 2fc ) = ^Afe - a<S> 2k , fc > 1. 



Therefore, 



-\/\n) - a$n-i, neven 
\JX( n ) - a$ n+ i, nodd 



= (-l) rt+1 y / ^(n> -a $„_(_!)«, nGN, 

with the convention that $_i = 0. To finish the proof it is now sufficient to observe that 
a double application of D maps, up to a multiplicative constant, $ n onto itself, 

D 2 $ n = — (L - a)$ n = -(A(„) - a)$„. 

□ 

Note that here, in contrast with the examples considered in [3] Section 7], Dj has 
proper invariant subspaces that decompose orthogonally the whole subspace U L 2 (X, /x) C 
L 2 (X, fi). They are spanned by the pairs {$„, $ n _ ej }, where n is such that rij > is even. 
Notice that acts trivially on the subspace spanned by {$ n : rij = 0}. 

Corollary 4.3. Gzwen / G N d \{(0, . . . , 0)} we have 

D%i = ( _ 1) H[/*f|(^/ 2 )ii( r JJ ( A W - a,)^ V M) „ r , n G K 

where I is a multi-index such that lj = if lj is even and lj = 1 otherwise, (— l) n l = 
{{-l) ni h, (-l)n*) and (n + 3/2)7= ((m + 3/2)^, . . . , {n d + 3/2)I d ). 

As a consequence of the above corollary and Bessel's inequality we get the following. 

Proposition 4.4. The series defining the Riesz transforms R converge in L 2 (X, fi) and 
for each order N > 1 the mapping 



\l\=N 



1/2 

2 



zs a (nonlinear) contraction in L 2 (X, /i). In particular, each R is a linear contraction. 

It is remarkable that the present approach to the higher order Riesz transforms is consid- 
erably simpler than that in [3J Section 4]. This is due to the fact that in the symmetrized 
setting the subspace spanned by the orthogonal system is invariant under actions of the 
associated 'derivatives'. More comments in this connection will be given in Section [51 

We pass to defining conjugate Poisson integrals in the symmetrized setting. The Poisson 
semigroup {Pt}t>o associated with £ is, by the spectral theorem, given on L 2 (X,/i) by 

oo 

PJ = exp ( - t^ 2 )f = exp(-tA^ 2 )7W. 

m=0 
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Clearly, each P t , t > 0, is a contraction on L 2 (X, fi). We now define the conjugate Poisson 
integrals Uf, t > 0, j = 1, . . . , d, as the contractions on L 2 (X, /i) given by 

Uif = P t R j f, feL\X,fi). 
To rewrite this by means of the spectral series observe that by Lemma 14.21 



R jf = E ( (n) i : - ) (/• $ «><^ 



(n) 



M^n— (-1) iej 



Since (n — (— f)'™ 3 ej) = (n), we see that 

/ x /a w -a \ 1/2 

^/ = E ("l) ni+1 eX P -tv/A^ i (/,<fn)^n-(-ip er 

This, together with Bessel's inequality, shows that for each t > also the mapping 



/-> viwi 2 + ... + iwi 2 

is a contraction in L 2 (X, /x). 

Our definitions of Riesz transforms and conjugate Poisson integrals are well motivated 
by the following system of Cauchy-Riemann type equations. 

Proposition 4.5. Let f belong to the subspace of L 2 (X, //) spanned by the <& n 's. Then 

D t Uif = D 3 U l J, i,j = l,...,d, 

D 3 Ptf = -^U?f, j = l,...,d. 



If A = 0, then also 



for A > the function f on the right-hand side above must be replaced by f — A£ l Hof ■ 
Moreover, we have the harmonicity relations 



J^-L)iV = 0, (|L-L)C£/ = 0, j = l,...,d. 



t . 

d 



Proof. It is enough to restrict the situation to / = $ n , n e N . Since 

A (i) -a.X 1 / 2 

f'7'I\, = l-ir +1 exp (-^xAm) ( ^ 3 ) 



and for z 7^ j, A$n-(-i)"j ej = (- 1 ) n<+1 ( A (l ) > j - ai) 1/2 $n-(-i)"< ei -(-i) n i ej , the first identity 
follows. The second identity may be also easily justified because 

D^ n = (-ip +1 exp (-tyfi^ (A^ - a,) 1/2 $ n _ ( _ ipej = ~TJ{^ n . 

To verify the third identity, we observe that 

d 



Q^Pt^n = - V A(n> exp ( -£^/ A( n ) ) $ r 
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and since Dj commutes with L, thus also with Ul, we have 

DjVi* n = ui(D^ n ) = (-i)»i+yA (n) . - aj ui$ n _ { _ irjG . 

= ~{ X (X - ^ 6XP (- t \A^) {\n)Y 1/2 <bn- 

Here the last equality is obtained by recalling that (n — (— l) nj e,,-) = (n) and also noticing 
that n — (—l) nj ej — (— 3 ej = n. Finally, checking the harmonicity relations does 
not cause any problems. □ 

We remark that a suitable information on the growth of the eigenvalues A/ n \ and on the 
growth of the eigenfunctions $ n and their derivatives allows to show that the identities of 
Proposition 14.51 hold in fact for all / G L 2 (X, /i); see [3j Proposition 5]. 

Further support for the symmetrized conjugacy scheme is provided by the identity 

d 

Y,R)f = -/ + AC" 1 /, / e n L 2 (x,u) ; 

notice that when A = the potential term above vanishes. This is an analogue of the well- 
known relation J^ . R? = —Id, satisfied by the classical Riesz transforms Rj = <9j(— A) -1 / 2 . 

A comment concerning the 'complex' picture from Section [3] is in order. Note that 
replacing the symbols N d , (n) and $ n in (14.11) and (14. 2 j) by Z d , |n| and ^/ n , respectively, 
changes neither Dom£ nor £. Further, replacing A( n ) and $ n in (14. 3 j) by A| n | and ^/ n , 
respectively, does not change the Riesz operators (in the one-dimensional setting, the action 
of £) on \l/ n is D^ n = isgnn^Aini — a^ n , n G Z, and similarly for -D^). Consequently, 
the Poisson semigroup and the conjugate Poisson integrals remain unchanged. 

Finally, notice that in the context of Example 1 the Riesz transform given by (14. 3 j) for 
I = d = 1 results in the classic conjugacy mapping 

oo oo 

— + (a n cos nx + & n sin nx) h- >■ (6 n cos nx — a n sin nx) 

n=l n=l 

(Enez a n emi ^ Snez i s S n n a n emx 111 * ne 'complex' picture). 

5. Comments and examples 

First we observe that the setting considered in [3] is naturally embedded in the sym- 
metrized situation. Indeed, given a function / on X, consider its extension / to X that 
is even with respect to all coordinates. Then the definitions and relations from the sym- 
metrized scheme can be applied to /, and this clearly induces analogous restricted defini- 
tions and relations related to the original space X . In this way the general definitions of 
Riesz transforms of order one given in [31 Section 3] and conjugate Poisson integrals given 
in [3], Section 5] are contained in the symmetrized definitions from Section |U In a similar 
manner the Cauchy-Riemann type equations and harmonicity relations [31 (5.3)-(5.6)] are 
'projections' of the identities from Proposition 14.51 Moreover, by considering extensions 
of a function / that are odd with respect to one coordinate and even with respect to all 
remaining coordinates it can be seen that the 'supplementary' operators and relations es- 
tablished in [21 Section 6] are suitable 'projections' of the symmetrized counterparts from 
Section HI 
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However, the definition of higher order Riesz transforms induced by the symmetrized 
scheme in the initial setting is essentially different from that postulated in [21 Section 
4]. Nevertheless, it seems to be far more appropriate and natural. Observe, that the 
'projection' from the symmetrized situation via considering functions that are even with 
respect to all coordinates leads to higher order derivatives in the initial setting that are of 
the form 



where n = (m, . . . , n d ) is a multi- index. This obviously makes a contrast (when rij > 1 for 
some j = 1, . . . , d) with the derivatives 



used in [3] to define higher order Riesz transforms. 

The definition of higher order Riesz transforms in the initial setting based on D n (i.e. 
the 'even projection' of the symmetrized definition) seems to be more natural, in particular 
no complications occur in connection with showing L 2 -boundedness of these operators, see 
[21 Section 4, Section 7.9]. The new light on understanding higher order derivatives and 
Riesz transforms in the initial setting should also have an important impact on developing 
the theory of Sobolev spaces related to orthogonal expansions. This subject remains to be 
investigated. 

We conclude the paper with several concrete examples involving selected classical or- 
thogonal expansions, where the symmetrization procedure can be easily traced explicitly. 
More exemplifications can be derived from those given in j3j Section 7]; in particular, we 
follow the notation from there. For the sake of clarity and simplicity, in Examples 2-4 
below we assume that d — 1. 

Example 2. Let {h n : n e N} be the classical Hermite functions on R and consider the 
system = V / 2^2fc on the half-line X = (0, oo), iel This system is an orthonormal 
basis in L 2 (X,dx) consisting of eigenfunctions of the harmonic oscillator L = —-jL + x 2 
restricted to (0, oo). The related derivatives decomposing L are 6 = -^+x and 6* = 
see [21 Section 7.4] (notice that 6 is not skew-symmetric). Passing to the symmetrized 
situation we get the orthonormal system {$ n } in L 2 (R,dx), which coincides, up to signs, 
with the full system of Hermite functions. The symmetrized derivative is Df = & + xf, 
where f(x) = f(—x) is the reflection of /, and the symmetrized 'Laplacian' has the form 



oscillator by the reflection term above. On the other hand, the derivative D is formally 
skew-adjoint in L 2 (IR, dx). 

Example 3. A natural generalization of the previous example is obtained by taking 
X = (0, oo) equipped with the measure fi a (dx) = x 2a+1 dx, a > —1, and considering the 
system ip k = of Laguerre functions of convolution type, see [21 Section 7.6]. Here a is 
a parameter of type, and the value a = —1/2 corresponds to the situation described in 



©"=(... 5,616,616, )(... 6 2 5* 2 5 2 6* 2 5 2 )...(... 6 d 6* d 6 d 6* d 6 d ) , 
\ s ,i \ v ,i \ s ,/ 




ni components ni components rid components 



6 n = 6^ ... 5 n d 



with / odd 
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Example 2. The related standard 'Laplacian' is 



d 2 2a + 1 d o 
+ x l 



dx 2 x dx 

and the associated derivatives are of the form 6 = 4- + x, 5* = + 2 ^- s ^. Passing to 
the symmetrized situation we arrive at the system {$ n } that coincides, up to signs, with 
the system of generalized Hermite functions emerging in the context of the Dunkl harmonic 
oscillator and the underlying group of reflections isomorphic to Z2, see jl]. However, the 
symmetrized 'Laplacian' 

d 2 f 2a + l df 2 2a + l 

= —t-= -r- + X f + — /odd + 2/ dd 

ax z x dx x 1 

differs from the Dunkl harmonic oscillator by the term 2/ D dd above. The symmetrized 
derivative Df = ^ + xf + ^^-/odd is skew-symmetric, which is not the case of 5. 

Example 4. Finally, consider an orthonormal basis {(fk} of L 2 (X,n) consisting of eigen- 
functions of a divergence form operator 



t f 1 f f\' 1 t d2 f w ' d f , t 
Lf = {wf ) +af = --— T + *f, 

w dx z w dx 

where w is the density of /1 and a > is a constant. We assume that all the technical 
assumptions from Section [2] are satisfied, in particular X is an interval of the form (0,c), 
< c < 00. The derivatives decomposing L have the form 5 — 4-, 5* = — -r- — — . 
Performing the symmetrization procedure, we find the symmetrized 'Laplacian' 

w d *f w ' d f ^ f ( w '\'t 
dx z w dx \w J 



Df = + + -/odd, 



and the associated derivative 

dx w 

which is skew-symmetric in L 2 (X, /i). 

A special case of the situation just described occurs when tp k are the (normalized) 
Hermite polynomials of successive even orders, X = (0, 00), w(x) = e~ x , a = 0, and 
L = —4^2 + is the classical Ornstein-Uhlenbeck operator restricted to the positive 
half-line. Passing to the symmetrized situation one receives the (normalized) system of 
Hermite polynomials of all successive orders and the symmetrized 'Laplacian' L which 
differs from the Ornstein-Uhlenbeck operator by the reflection term 2/ Q dd- The point, 
however, is that the associated derivative is skew-symmetric. 

Another important special case is obtained by choosing ip^ to be the normalized Jacobi 
trigonometric polynomials considered on the interval X = (0, tt) equipped with the measure 
dfj,(9) = w(9)d9 = (sin |) 2a+1 (cos ^) 2/3+1 d9. Here a, f3 > — 1 are parameters of type, and 
taking a = (3 = — 1/2 we recover the situation of cosine expansions already discussed in 
Example 1. The related 'Laplacian' is 

d 2 9 a- [3 + (a + (3 + l)cos9 d ^a + /3 + l\ 2 
~ ~d9^ sin~0 d9 + 
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and the derivatives decomposing it have the form 5 — 5* — — — (a + \) cot § + (/? + 
|tan~). The symmetrized 'Laplacian' is then 

d 2 f a-f3 + (a + (3 + l)cos6 df / a + (3 + i y 

J ~ d9 2 sm9 d9 + \ 2 ) J 

(a + ft + 1) + (a - f$) cos6 

H • 2 a J odd 

sin v 

and the associated skew-symmetric derivative is 

df a- /3 + (a + (3 + 1) cos6 

u 3 — -77H J° dd - 

de sm y 

It is worth to note that this D coincides with the Jacobi-Dunkl operator on the interval 
(— 7r, 7r), and the extended system in the 'complex' picture {^ n } consists of trigonometric 
polynomials called the Jacobi-Dunkl polynomials; see [TJ, for instance. 

Developing widely understood harmonic analysis for orthogonal expansions is intimately 
connected with (sometimes implicit) choice of the associated 'Laplacian'. The results of 
this paper, and in particular the examples given above, show that in many cases there are 
reasonable and in some aspects more natural alternatives for standard 'Laplacians' related 
to various orthogonal systems appearing in the literature. From this point of view deleting 
the constant A in the decomposition (13. ip of L would lead, in some sense, to canonical 
'Laplacian' associated to general orthogonal expansions, bringing the related harmonic 
analysis closer to the classical case. 
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